Astronomy & Astrophysics manuscript no. corotexo2 


© ESO 2009 


September 22, 2009 





Properties of starspots on CoRoT-2* 

Adriana Silva-Valio\ A. F. Lanza^, R. Alonso^, and P. Baige'' 



o 



' CRAAM, Mackenzie University, Rua da Consolajao, 896, 01302-907, Sao Paulo, Brazil 

^\ ' e-mail: avalio@craam.mackenzie .br 

, ' INAF-Osservatorio Astrofisico di Catania, Via S. Sofia, 78, 95123 Catania, Italy 

■ ^ Laboratoire d'Astrophysique de Marseille, UMR 61 10, CNRS/Universite de Provence, Traverse du Siphon, 13376 Marseille, France 

(N ■ 

^ ' Received ; accepted 

; ABSTRACT 

00 : 

' Context. As a planet eclipses its parent star, a dark spot on the surface of the star may be occulted, causing a detectable variation in 

' the light curve. 

Aims. Study these light curve variations during transits and infer the physical characteristics of the stellar spots. 
I— I . Methods. A total of 77 consecutive transit light curves of CoRoT-2 were observed with a high temporal resolution of 32 s, correspond- 

' ing to an uninterrupted period of 134 days. By analyzing small intensity variations in the transit light curves, it was possible to detect 

' and characterize spots at fixed positions (latitude and longitude) on the surface of the star. The model used simulates planetary transits 

• , and enables the inclusion of spots on the stellar surface with dilferent sizes, intensities {i.e. temperatures), and positions. Fitting the 

. data by this model, it is possible to infer the spots physical characteristics. Because what is observed is the stellar flux blocked by the 

(Dlh' spots, there is a degeneracy between the spots intensity and area. Thus the fits were either in spot longitude and radius, with a fixed 

intensity, or in spots longitude and intensity, for spots of constant size. The model allowed up to 9 spots to be present at the stellar 
surface within the transit band. 

Results. Before the modeling of the spots were performed, the planetary radius relative to the star radius was estimated by fitting the 
^ ' deepest transit to minimize the elfect of spots. A slightly larger (3%) radius, 0.172 R„„r, resulted instead of the previously reported 

I I ' 0.1667 R„ar. The fitting of the transits yield spots, or spot groups, with sizes of ranging from 0.2 to 0.7 planet radius, Rp, with a mean 

of 0.41 ± 0.13Rp (~ 100,000 km), resulting in a stellar area covered by spots within the transit latitudes of 10-20%. The intensity 
. varied from 0.4 to 0.9 of the disk center intensity, 4, with a mean of 0.60 ± 0.19/,., which can be converted to temperature by as- 

^ ' suming blackbody emission for both the photosphere and the spots. Considering an elfective temperature of 5625 K for the stellar 

Ijr^ ' photosphere, the spots temperature ranges mainly from 3600 to 5000 K. 

l/^ , Conclusions. The spot model used here was able to estimate the physical characteristics of the spots on CoRoT-2, such as size and 

. intensity. These results are in agreement with those found for magnetic activity analysis from out of transit data of the same star. 

, Key words, extra-solar planets; star spots; stellar magnetic activity 
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1. Introduction The n iagnetic activity of CoRoT-2 star has been studied in 

r. detail by iLanza et"aLl (l2009l) who analyzed its out of transit 

>► Exacdy four centuries ago, spots on the surface of the Sun were ijght curve with modulations of -6% of its total flux. Using 

first detected by Galileo. Sunspots are cool regions of strong maximum entropy regularized models, Lanza and collaborators 

1.^ concentration of magnetic fields on the photosphere of the Sun. (2009) modeled the light curve considering both the presence of 

Presently, thanks to the observations with unprecedent e d pho - sunspots and faculae (dark and bright regions, respectively). This 

■ - - ' tometric precision of the CoRoT satellite paglin et alj |2006|), study detected the existence of two active longitudes located in 

spots on the surface of another star can be detected, that of opposite hemispheres, and also that these longitudes varied with 

CoRoT-2, a much more active star Such activity is proba- time. The active longitudes are regions where the spots prefer- 

bly a conseq uence of its young age, estimated in 0.5 Gyr ably form. The total area covered by the spots were seen to vary 

dBouchy et aU 2008). CoRoT-2 is one of the 7 stars with tran- periodically with a period of 29 + 4 days. The authors were also 

siting planets detected so far with CoRoT. able to estimate an ampUtude for the relative differential rotation 

Spot activity has been infeiTed in the past from the modu- of < 0.7%. 

lation observed in the light of stars. As the star rotates, differ- Here we propose to study these same spots, however, with 

ent spots on its surface are brought into view. Because spots are a different approach. A transiting planet can be used as a probe 

cooler, and therefore darker, than the surrounding photosphere, of the contrasting features in the photosphere of the star When 

the detected total light of the star diminishes as different spots the planet occults a dark spot on the surface of its ho st star 

of varying temperature and size face the observer This periodic small variations i r i the light curve ma y be detected (ISilval 120031; 

modulation enables the determination of the rotation period of |Pont et all 120071; ISilva- Valid |2008|) . From modeling of these 

the star. variations, the properties of the starspots can be determined, such 

as size, position, and temperature. Moreover, from the continu- 

* CoRoT is a space project operated by the French Space Agency, ous and long duration of observation provided by the CoRoT 

CNES, with participation of the Science Programme of ESA, satellite, the temporal evolution of individual spots can be ob- 



ESTEC/RSSD, Austria, Belgium, Brazil, Germany, and Spain. tained. 
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The next session describes the observation performed by 
CoRoT, whereas the model used here is introduced in the follow- 
ing session. The main results are presented in Session 4. Finally, 
a discussion of the main results and the conclusions are listed in 
Session 5. 

2. Observation of CoRoT-2 

A planet around the star CoRoT-2 was detected during one of 
the long run observations of a field toward the Galactic cen- 
ter performed by the CoRoT satellite, and is the second transit- 
ing planet discovered. The planet, a hot Jupiter, with a mass of 
3.3 Mjup and radius of lA lRjuo, orbits its host star in just 1.73 
day (lAlonso et al.L l2008l) . CoRoT-2 is a solar- like star of type 
G7, w ith 0.97 Mp and 0.90 2 Rq that rotates with a period of 4.54 
days (lAlonso et al.L |2008l) . The parameters of the star and the 
planet plus the orbital ones were deter mined from a com bina- 
tion of CoRoT photometric light curve (lAlonso et al.L l2008i) and 
ground based spectroscopy (iBouchv et al.Ll2008lK 

The data analyzed here were reduced following the same 
procedure as Alonso et al. (2008|). The light curve was filtered 
from cosmic ray impacts and orbital residuals. Then this cleaned 
light curve was folded considering an orbital period of 1.743 
day, from which the parameters of the planetary system were 
derived. Beside the planet orbital period of P,, = 1 .743 day and 
stellar rotation of Pj - 4.54 days, the orbital parameters were: 
semi-major axis of a = 6.7 star radius (Z?,) and inchnation angle 
of 87.84°. 

A total of 77 transits were detected in the light curve with a 
high temporal resolution of 32 s, in a total of 134 days. The rms 
of this signal was estimated from the out of transit data points to 
be cr = 6 X 10 '* in relative flux units. Small light variations were 
detected during the transits, usually with fluxes between 3 and 10 
sigma above the transit model without spots, the largest variation 
reaching 18 sigma. These intensity variations were interpreted 
as the signatures of star spots and were thus modeled. Figure [T] 
shows the light curves from all transits, where the vertical extent 
of each data point represents the rms of the signal. Also plotted 
on this figure is the model transit considering that no spots are 
present on the stellar surface (gray curves). 

3. The model 

The physical characte ristics of star spots is obtained by fitting 
the model described in lSilval ( l2003h . In the model adopted here, 
the star is a 2-D image with intensity de creasing following a 
quadratic limb darkening law according to lAlonso et al.l (l2008l) . 
whereas the planet is taken as a dark disc. The modeled light 
curve of the transit is obtained as follows. The planet position 
in its orbit is calculated every two minutes and the total flux is 
just the sum of all the pixels in the image (star plus dark planet). 
This yields the light curve, that is, the relative intensity as a func- 
tion of time during the transit. The model assumes that the or- 
bit is circular, that is, null eccentricity (consistent with the mea- 
sured eccentricity of 0.003 + 0.003), and that the orbital plane is 
aligned with the star equator. In the case of CoRoT-2 the latter is 
a good assumption since, by measuring the Rossiter-McLaughlin 
effect, iBouchv et alj (l2008b obtained that the angle between the 
stellar rotation axis and the normal of the orbital plane is only 
7.2 + 4 .5°. The orbital parameters were taken from lAlonso et al.l 
( |2008|) . such as period of 1.743 day and orbital radius of 6.7 

The model also allows for the star to have features on its 
surface such as spots. The round spots are modeled by three pa- 




Fig. 1. All the 77 light curves during the transit of CoRoT-2b 
in front of its spotted host star The gray solid line represent the 
model of a spotless star. The error bars of the data points indicate 
therms of 0.0006. 

rameters: (i) intensity, as a function of stellar intensity at disk 
center, 7^. (maximum value); (ii) size, or radius, as a function 
of planet radius, Rp\ and (iii) position: latitude (restricted to the 
transit path) and longitude. On all the fitting runs, the latitude 
of the spots has remained fixed and equal to the planetary tran- 
sit line, which for an inclination angle of 87.84° is -14.6°. This 
latitude was arbitrarily chosen to be South, thus the minus sign. 
The longitude of the spot is measured with respect to the central 
meridian of the star, that coincides with the line-of-sight and the 
planet projection at transit center. 

When the spot is near the limb, the eff'ect of foreshortening is 
also featured in the model. However, this model does not account 
for faculae. Solar-like faculae have a negligible contrast close 
to the disc center, having significant contrast only close to the 
limb. Si nce the facular-to-s potted area ratio for CoRoT-2 is only 
Q=1.5 jLanza et al 1 120091) instead of 9 as in the Sun, and we 
limit our analysis to spots between -70 and -1-70" from the central 
meridian, the photometric eff'ect of faculae is very small and can 
be safely neglected. 

The blocking of a spot by the planet during a transit implies 
an increase in the light detected, because a region darker than the 
stellar photosphere is being occulted. Thus the effect of many 
spots on the stellar surface is to decrease the depth of the transit, 
as can be seen in Figure [T] This is turn will influence the de- 
termination of the planet radius, since a shallower transit depth 
results in a smaller estimate o f the planet diameter when spots 
are ignored ( Silva-VahdlloTol) . 

In this work, to estimate the best model parameters for the 
planet and its orbit, the deepest transit was sought. Of all the 77 
transits, the 32nd transit displays the lower light variation (see 
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Fig. 2. Top: Examples of transit light curves of CoRoT-2. The 
black crosses represent the 32nd transit assumed to occur when 
the star has minimum spot activity within the transit latitudes, 
whereas the gray data points represent a typical transit (the fifth 
one). The model light curve without any spots is shown as the 
thick solid curve. Bottom: Residual from the subtraction of the 
data minus the model for a star with no spots for the 5th (gray) 
and 32nd (black) transits. The vertical bars in both panels repre- 
sent the estimated uncertainty of 0.0006 from out of transit data. 



star f' 


«) ' - 

--^ spot 




1 

I planet,,,,--''''''^ 







-3 3 

. (».„..) 




0.96 E ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ : 

-4-2 2 4 

lim.(h) 

Fig. 3. Top: Top view of the star and the planet in its orbit. A 
spot at 45° longitude on the stellar surface is depicted. Bottom: 
Light curve obtained from a star with a spot at 45" longitude. The 
dotted line represents a model for a transit in front of a spotless 
star 



Figure [TJ. This was interpreted as the star having the minimum 
number of spots on its surface within the transit latitudes during 
the whole period of observation (134 days). The 32nd transit 
is shown in Figure |2] as black crosses, for comparison, the fifth 
transit (gray crosses) is also shown in the same figure. 

A light curve obtained from a model without any spot is 
shown as a thick solid curve on the figure. However, to obtain 
this model light curve it was necessary to use a planet radius of 
0.172 stellar radius, instead of th e 0.1667 stellar radius quoted 
on Table 1 of lAlonso et al. I (I2008h . an increase of about 3%. We 
note that this may not be a difference in the actual radius but 
rather an artifact due to the uneven spot coverage on the total 
surface of the star during that specific transit. Nevertheless, this 
implies that star spots can hinder the exact size estimate of a 
planet by ma king the transit hgh t curve shallower than it would 
otherwise be ( ISilva-Valioll201(]h . 

The actual radius of the planet may very well be 0. 1667 Rs,ar, 
since this was calculated from phase folded and averaged light 
curve. Supposing that the average area of the star covered by 
spots does not change during the whole period of observation 
(134 days), then when there are few spots along the transit line 
band {e.g. 32nd transit), there should be more spots on the re- 
mainder of the star 

3.7. Spot modeling 

As mentioned in the previous section, the spots can be modeled 
by three basic parameters: intensity, radius, and longitude (the 
latitude is fixed at -14.6°). All the fit s were performed using the 
AMOEBA routine (iPress et al.Lll992l) . The longitude of the spot, 
is defined by the timing of the light variation within the transit. 
For example, the small "bump" seen in the fifth transit (gray 
crosses in Figure |2]i slightly to the left of the transit center, at 
approximately 0. 1 h, is interpreted as being due to the presence 
of a spot at a longitude of 5.6°, where longitude corresponds 



to the line-of-sight direction, taken as the central meridian of the 
star According to the diagram shown in Figure [3j the longitude 
of a spot may be estimated as: 

/ a \ (f/24) 
6 - sm smy6 — cos a where /3 - zn (1) 

\ I Po 

where t is the time, measured with respect to the transit center, 
and given in hours, P„ is the orbital period in days, and a is the 
latitude of the transit. The above equation is used to estimate 
the initial guess of the spots longitude, which was one of the 
parameters to be determined from the model fit to the data. 

The next step was to decide the maximum number of spots 
that were needed to fit each transit. Models with a maximum 
of 7, 8, and 9 spots on the stellar surface during each transit 
were tried. The results obtained from each run were qualitatively 
similar in all cases. It was found that in the case of 9 spots, the 
residuals were smaller than the uncertainty of the data (0.0006). 
Therefore, the results reported here are those from the fits with 
up to a total of nine different spots per transit. 

4. Spot parameters resulting from the fits 

CoRoT-2 is a very active star, and many intensity variations were 
identified in each transit, implying that there are many spots 
present on the surface of the star at any given time. Beside its 
longitude, the spot signature also depends on its intensity and 
size. In fact, the flux perturbation from the spot is the product 
of the spot intensity and its area. Thus there is a degeneracy be- 
tween the values of the radius and the intensity of one spot. 

Thus, to minimize the number of free parameters in the mod- 
eling, there are two alternatives. The first approach is to fix the 
radius of all the spots, for example, as half the planet radius, and 
fit each spot intensity and longitude for all transits. Another way 
is to fix the intensity of all the spots with respect to the maximum 
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Fig. 4. Top: Two examples of the model of the synthesized star 
with spots for transit 13. The left panel shows spots with fixed 
radius of 0.5 Rp and varying intensity, whereas the right panel 
presents spots with constant intensity of 0.5 Ic and changing ra- 
dius. Bottom: Residuals of the two models of constant radius 
(black) and intensity (red). 



central intensity, 4, and allow the spot radius and longitude to 
vary in each fit of the transit light curve. 

Examples of the fitting by the two methods are shown in 
FigurelH The two top panels represent the synthesized star with 
spots of: varying intensity and fixed radius of 0.5 Rp (left) and 
varying radius and fixed intensity at 0.5 /(■ (right). The residuals 
of the data minus the two fits are shown in the bottom panel and 
show the little difference between the methods. 

The star was considered to have a varying number of spots 
(from 2 to a maximum of 9), with spot position defined at a cer- 
tain longitude but a constant latitude of -14.6° (the transit lati- 
tude). The longitude considered is the topocentric longitude, that 
is, zero angle is defined as that of the line-of-sight, or transit cen- 
ter, when star, planet, and Earth are aligned. 

Figure |5] shows the total relative flux calculation, that is, the 
sum of all spot contrast times the squared radius (or area) in 
each transit, for both models. The spot contrast is taken as 1 - f,, 
where is the relative intensity of the spot with respect to disc 
center intensity. As mentioned above, the flux, F, of a single spot 
is: F oc (l - fi)R^pgf. For each transit the total relative flux was 
calculated by summing the flux of individual spots. As expected, 
the results from both models agree. 

4.1. Spot longitudes 

A histogram of the spots longitudes is shown in Figure|6]for the 
two models considered. Basically, seven main longitudes may be 
identified in the figure. Fortunately, they are approximately the 
same seven longitudes on both models, which is reassuring. The 
figure shows a slight predominance of the spots location at zero 
longitude, that is, the angle in the direction of the planet. 

These are topocentric longitudes, that is, they are not the 
ones located on the rotating frame of the star, but rather are 
measured with respect to an external reference frame. In or- 
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Fig. 5. Total relative flux of all spots per transit for both the 
model with constant spot intensity, 0.5 4 (black), and constant 
radius, 0.5 Rp (gray). 
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Fig. 6. Histogram of spot longitudes for both models with con- 
stant intensity (top) and radius (bottom). The vertical dotted lines 
represent longitudes of -67, -43, -18, -1, 25, 42, and 67°. 



der to obtain the longitudes in the st ellar rotating f rame, one 
needs an accurate period for the star 'Alonso et al.' (20081) re- 
port a 4.54 days period, whereas Lanza et al. ( 2009) obtained a 
period of 4.52 + 0.14 days in order to fit the rotational mod- 
ulation of the out of transit data. A more precise estimate of 
the period is need in order to analyze the spot results. A de- 
tailed investigation of the rotational longitudes and the spot life- 
time is und erway and will be reported in an accompanying paper 
(ISilva-Valio & LanzaLl20ld) . 

4.2. Spot intensity and temperature 

Spots with smaller intensity values, or high contrast spots, are 
spots cooler than those with intensity values close to Ic. The spot 
intensities obtained from the model with spots of fixed radius of 
Q.SRp are shown in the top panel of Figure |7] The figure shows 
that the spot intensities range from 0.4 to 0.9 Ic, with an average 
value of 0.60 + 0.19 of the stellar central intensity. This value is 
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Fig. 7. Results from the fit of the model with spots of fixed radius 
for all transits. Distributions of the spot intensity as a fraction of 
the central intensity (top), and spot temperature (bottom). 



close to the value of 0.665/^ used bv lLanza et al.l (l2009l) . that is 
the mean value of spots on the Sun. 

These intensities can be converted to spot temperature by 
assuming blackbody emission for both the photosphere and the 
spots. The temperature is estimated as: 



spot 



hv 
k'b 



In 



1 + 



expl'^fV// -1 



fi 



where Kb and h are Boltzmann and Planck constants, respec- 
tively, V is the frequency associated to a wavelength of 600 nm, 
fi is the fraction of spot intensity with respect to the central stel- 
lar intensity, 4, and T^ff is the effective temperature of the star 
Considering T^f j - 5625 K ( Alonso et al., 2008), the spots tem- 
perature range from 3600 to 5500 K, which are 100-2000 K 
cooler than the rest of the disk. The mean temperature of con- 
stant size spots on CoRoT-2 is 4600 + 400 K. 

4.3. Spot radius 

The distribution of spot radius from all the transits resulting from 
three simulations with different spot constant intensities: 0.3, 
0.5, and 0.665 of the central intensity, Ic are shown on Figure|8] 
As can be seen from the histograms, as the spot constant inten- 
sity increases (or conversely its contrast decreases), the resulting 
radius also increases. This occurs in order to keep the spot flux 
the same. 

The models with fixed spot intensity at 0.3, 0.5, and 0.665/f 
resulted in spots with average radius of 0.34 + 0.10 Rp, 0.41 + 
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Fig. 8. Results from the model fit for all transits: distributions of 
the spot radius in units of planet radius for three different fixed 
intensities as a fraction of the central intensity: 0.3 4 (top), 0.5 
Ic (middle), and 0.665 /e(bottom). 



0.13 Rp, and 0.50 + 0.17 R,,, respectively. The 0.5 R,, value as- 
sumed in the spot model with fixed radius agrees with the aver- 
age radius of spots with 0.6654, w hich interestingly is the mean 
intensity of sunspots lLanza et al] (i2009.) . 

The results show that the radius of the modeled spots varies 
(-2) from 0.2 to 0.7 Rp. Assuming a planet of 1 .4657? y„p, this implies 
in spots with diameters of 40 to 150 Mm, with a mean value of 
about 100 Mm. 



4.4. Stellar surface area covered by spots 

To estimate the area covered by spots on the surface of the star, 
only the significant spots should be taken into account. By sig- 
nificant we mean spot with contrast larger than 10% of the stel- 
lar central intensity or radius larger than O.lRp, depending on 
the model. According to this criterion, the number of spots on 
the surface of the star during each transit varied from 2 to 9 
spots, with an average of 7 (for the constant radius model) or 
8 (model with constant spot contrast) spots per transit. The tran- 
sit with smaller number of spots corresponds to the 32nd transit 
discussed before. Histograms of the number of spots per transit 
are shown in the two top panels of Figure|9]for models with spots 
of constant radius (left) and intensity (right). 

The mean surface area covered by the spots during each tran- 
sit is the sum of the area of all spots detected in that transit di- 
vided by the total area occulted by the passage of the planet (see 
the diagram on Figure[TO|. The total area, A/gi, of the transit band 
occulted by the planet is calculated as: 
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Fig. 9. Top: Distribution of the number of spots detected on each 
planetary transit for both the constant radius, 0.5 Rp (left), and 
intensity 0.5 Ic (right) spot models. Bottom: Histograms of the 
stellar surface area covered by spots within the transit latitudes 
for each model (constant radius and intensity, left and right pan- 
els, respectively). 
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Fig. 10. Diagram of the total area of the star occulted by the 
planet transit, represented by the hatched region. 




Fig. 11. Mean stellar surface area covered by spots within the 
transit latitudes as a function of fixed spot intensity. 
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The 7/9 factor arises because in the fit, only the spots between 
longitudes -70° and +70° are considered due to the difficulties 
in fitting spots too close to the limb where the light curve is very 
steep. 

The surface area covered by spots on each transit, taking 
into account only the area of the transit band, of course, com- 
puted from the above equation is show in the bottom panels of 
Figure |9] for both models. For the model with constant radius 
(bottom left), the spot surface area for each transit is the product 
of the number of spots on that transit multiplied by ji{Q.5Rp)^. 
On the other hand, for the model with constant spot intensity, 
the star surface area covered by spots is the sum of the area of 
all spots, that is, E;r/?^^,^,, where Rspot is the radius of the spot 
obtained from the fits. The average values of the stellar surface 
area covered by spots are 16 + 3% and 13 ± 5% for the models 
with constant radius {Q.5Rp) and intensity (0.5/c), respectively. 

Because the data is only sensitive to the flux decrease due 
to the presence of spots, different values of spots intensity will 
produce fits with spots of different radius. The stellar area cov- 
ered by spots was also calculated for the models with different 
intensities (0.3, 0.5, and 0.665/c) discussed in the previous sub- 
section. A plot of the mean stellar area coverage as a function of 



the fixed spot intensity for each run is shown in Figure [TT] Also 
plotted on the figure as an asterisk is the area coverage obtained 
from the model with fixed radius at O.SRp, where the mean in- 
tensity of 0.60/c was considered. 

As can be seen from the figure, the area covered by spots 
within the transit latitudes increases as the spot intensity in- 
creases, varying, on average, between 10 and 20%. This occurs 
because to fit the same variation in flux, a hotter spot (less con- 
trast) needs to be larger, since the occulted area of the stellar 
surface is not so dark in this case. In summary, to account for 
the total decrease in Ught of the star due to the presence of hotter 
spots, one needs larger spots. This same trend was seen in the 
results of .Wolter et al.. C2009.) from modelling of a single spot. 



5. Discussion and conclusions 

CoRoT-2 star is a young and reasonably active star. The presence 
of spots on the surface of the star can influence the determination 
of the orbital parameters of a pl anet in orb i t by d istorting the 
transit light curve in two ways ( Silva- Valiol l2010h . One is the 
presence of spots on the limb of the star which will cause the 
transit duration to be shorter than it really is. The other distortion 
is to make the transit shallower if there are many spots on the 
surface of the star This would cause the planet radius estimate 
to be smaller than its real value. The latter eff'ect was observed 
in the dataset analyzed here, where the spot model applied here 
yi eld a radius of 0.17 2 R^tar instead of the 0.1667 Rstar listed 
in lAlonso et al.l (|2008'). In this case, we do not believe that this 
represents a real difference in the planet radius, but rather an 
artifact of the spot distribution on the surface of the star at a 
given time. 
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Here the star was modeled as having up to 9 round spots 
at any given time on its surface at fixed positions (latitude and 
longitude) during the uninterrupted 134 days of observation by 
the CoRoT satellite. For each transit, the longitudes of the spots 
were obtained from a fit of the model to the data. The other free 
parameter obtained from the fit was either the spot radius or its 
intensity. 

Two fitting approaches were performed on the transit light 
curves of CoRoT-2: the first one considered the radius of the 
spots to be fixed at 0.5 Rp, the second one kept the spot intensity 
at a constant value of the stellar central intensity, 7^. In the sec- 
ond approach, three different values of the fixed spot intensity 
w ere considered, one of them being the value of 0.665 Ic used 
in iLanza et alj (l2009l) . On every transit there were, on average, 
7-8 spots on the visible stellar hemisphere within longitudes of 
+70°. 

Despite the two methods with different fixed parameters of 
the spots, either radius or intensity, the results obtained for the 
spots characteristics were very similar. For example, the longi- 
tudes are approximately the same and the spot surface area cov- 
erage are compatible. Also, the model with spots of fixed inten- 
sity at 0.665/c yields a mean value for the radius of 0.5 1 ±0. 1 %Rp, 
agrees with the mean intensity of 0.60 + 0. 19/c obtained from the 
model with spots of constant radius, .5Rp. This mean in tensity 
is close to the value of 0.665 Ic used bv lLanza et al. (2009 ) which 
is the same as the sunspot bolometric contrast ( Chap man et al.L 
|1994|) . This agreement between the various approaches shows 
the robustness of the model. 

From the method of spots with fixed radius but varying in- 
tensities, the mean temperature of the spots was estimated as 
4600 ± 400 K by considering blackbody emission for both the 
stellar photosphere and spots. These spots are about 1000 K 
cooler than the surrounding photosphere. On the other hand the 
runs of spots with varying radius but fixed intensity of 0.3, 0.5, 
and 0.665 h yield mean radius of 0.35, 0.41, and 0.50 Rp. The 
increase in radius size means that darker spots (smaller intensity) 
are smaller than brighte r ones (intensities close to /c). 

IWolter et al.1 12009) modeled a single spot on one transit 
(here transit 54 of Figure [1]). The authors modeled the spot with 
different intensities, analogous to our procedure, and obtained a 
spot radius of 4.8° (in degrees of stellar surface) for 0.3/c. This 
specific feature, a "bump" on the light curve transit was better 
fit by our model by two spots at longitudes of 3 and 18°, with 
radius of 0.43 and Q39Rp for the model of constant spot inten- 
sity of 0.3/c. These radii are e quivalent to 4.2 and 3.9°, very 
similar to the results obtained bv l Wolter et alj (l2009l) . Moreover, 
IWolter et alj (2009) also confirm that spots with smaller contrast 
(higher intensity relative to the photosphere) need to be larger 

The spots on CoRoT-2, of the order of ~ 100,000 km, are 
much larger than sunspots, about 10 times the size of a large 
sunspot (10 Mm). The mean surface area of the star covered by 
spots within the transit latitudes is about 10-20%. This is larger 
than th e 7 to 9% of the total spotted area found bv [Lanza et alj 
(|2009|) . However, these values were estimated considering the 
whole star, that is, also the pole areas, where there are no spots 
in the case of the Sun. The values obtained here are only for the 
transit latitudes, which span approximately 20° and are close to 
the equator. In this case, the latitudes coincide with he so called 
royal latitudes of the Sun, where most of the sunspots occurs. It 
was reassuring to see that the results from both methods agreed 
very well with each other and specially with those of the out of 
transit data analysis Lanza et al. (2009). 

Long term observations such as the one provided by CoRoT 
are paramount to understand the physics of stars spots. The 



model applied here to the CoRoT-2 data is capable of obtain- 
ing the physical properties of spots with the advantage of fol- 
lowing thekjemjjoralero^ that is done in an accompanying 
paper dSilva-Valio & Lanzal 120101) . It will be very mterestmg to 
perform similar analyzes on data from other stars with plane- 
tary transits observed by the CoRoT satellite, especially for stars 
which are not solar-like. 

Acknowledgements. We would like to thank everyone involved in the planning 
and operation of the CoRoT satellite which made these observations possi- 
ble. A.S.V. acknowledges partial financial support from the Brazilian agency 
FAPESP (grant number 2006/50654-3). 



References 

Alonso, R., Auvergne, A., Baglin, A. et al. 2008, A&A, 482, L21 
Baglin, A. et al. 2006, 36th COSPAR Scientific Assembly, 36, 3749 
Blanter, E. M., Mouel, J.-L. Le, Perrier, F. & Shnirman, M. G. 2006, Solar Phys., 
237, 329 

Bouchy, P., Queloz, D., Deleuil, M., Loeillet, B., Hatzes, A. P et al. 2008, A&A, 
482, L25 

Chapman, G. A., Cookson, A. M., & Dobias, J. J. 1994, ApJ, 432, 403 

Hathaway, D. H. & Choudhaiy, D. R 2008, Solar Phys. 250, 269 

Lanza, A. E, Pagano, I., Leto, G., et al. 2009, A&A, 493, 193 

Petrovay, K. & van Driel-Gesztelyi, L. 1997, Solar Phys. 176, 249 

Pont, E, Gilliland, R. L., Moutou, C., Charbonneau, D., Bouchy, E, Brown, T. 

M., Mayor, M., Queloz, D., Santos, N. & Udry, S. 2007, A&A, 476, 1347 
Press, W. J., Teukolsky, S. A., Vetterling, W. T. & Elannery, B. P 

1992, Numerical recipes in EORTRAN, The art of scientific computing 

(Cambridge: University Press), 2nd. Edition 
Silva, A. V. R. 2003, ApJLetters, 585, L147 
Silva-Valio, A. 2008, ApJLetters, 683, L179 

Silva-Valio, A. 2010, in Proceedings of the lAU Symp. 264: "Solar and Stellar 
Variability Impact on Earth and Planets", 3-7 August 2009, Rio de Janeiro, 
Brazil (in preparation) 

Silva-Valio, A. & Lanza, A. E. 2010, A&A(in preparation) 

Wolter, U., Schmitt, J. H. M. M., Huber, K. E, Czesla, S., MuUer, H. M., 
Guenther, E. W. & Hatzes, A. P 2009, A&A, 504, 561 



